CLAIMS 

1. A non-invasive method of obtaining at least one 
predetermined target soundwave field in a substantially 
homogeneous medium (19) masked by a bone barrier (14) by 
causing sound signals to be emitted by at least one array 
(17) of transducers through said bone barrier (14) , 

the method being characterized in that it comprises: 
• a training stage comprising at least the following 
steps : 

la) making a three-dimensional image of the 
bone barrier (14) , at least in part by using X-rays, 
thereby obtaining a parameter representative of the 
porosity of said bone barrier at various points; 

lb) from said three-dimensional image, 
determining three-dimensional maps of at least density, 
soundwave speed, and soundwave absorption in said bone 
barrier; 

Ic) determining a specific position for the 
array (17) of transducers relative to the bone barrier 
(14) ; 

Id) simulating at least one propagation of 
soundwaves between at least one point of the 
substantially homogeneous medium and at least some of the 
transducers of the array of transducers on the basis of a 
mathematical model of propagation and said three- 
dimensional maps of density, soundwave speed, and 
soundwave absorption; and 

le) on the basis of said simulation, 
calculating individual sound signals to be emitted by at 
least some of the transducers of said array of 
transducers in order to obtain said target soundwave 
field; 

• and a stage of actually positioning the array of 
transducers on the brain barrier, this stage comprising 
the following steps: 
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2a) initially positioning the array (17) of 
transducers on the bone barrier (14) approximately in the 
specific position; 

2b) using at least some of the transducers of 
5 the array of transducers to perform echography to locate 
the position of said array (17) of transducers relative 
to the bone barrier (14) ; and 

2c) refining the position of the array (17) of 
transducers relative to said bone barrier as a function 
10 of the location determined in step 2b) . 

2. A method according to claim 1, in which, during step 
2b) , the relative position between said array (17) of 
transducers and the bone barrier (14) is located by 
15 determining an outside shape of at least a portion of the 
bone barrier by echography using at least a portion of 
said array of transducers, and comparing said outside 
shape with said three-dimensional image of the bone 
barrier . 

20 

3- A method according to claim 1, in which: 

• in step la) , an initial step is performed during 
which a locating device (8) is rigidly secured on said 
bone barrier (14), said locating device (8) being adapted 

2 5 to absorb X-rays; 

• during step la) , the three-dimensional image made 
also gives the position of the locating device (8) on the 
bone barrier (14) ; and 

• during step 2b) , the relative position between the 
30 array (17) of transducers and the locating device (8) is 

identified by echography. 

4, A method according to any preceding claim, in which 
the array of transducers is included in a fluid- filled 
35 tank (13) having at least one flexible wall (12) , and 
during the positioning stage, said flexible wall is 
pressed against the bone barrier (14) . 



5. A non-invasive method of obtaining at least one 
predetermined target soundwave field in a substantially 
homogeneous medium (19) comprising at least a portion of 
a brain masked by a bone barrier (14) comprising at least 
a portion of a skull surrounding said brain, by causing 
sound signals to be emitted by at least one array (17) of 
transducers through said bone barrier (14), 

the method being characterized in that it comprises: 
• a training stage comprising at least the following 
steps : 

la) making a three-dimensional image of the 
bone barrier (14) , at least in part by using X-rays, 
thereby obtaining a parameter representative of the 
porosity of said bone barrier at various points; 

lb) from said three-dimensional image, 
determining three-dimensional maps of at least density, 
soundwave speed, and soundwave absorption in said bone ' 
barrier; 

Ic) determining a specific position for the 
array (17) of transducers relative to the bone barrier 
(14); 

Id) simulating at least one propagation of 
soundwaves between at least one point of the 
substantially homogeneous medium and at least some of the 
transducers of the array of transducers on the basis of a 
mathematical model of propagation and said three- 
dimensional maps of density, soundwave speed, and 
soundwave absorption; and 

le) on the basis of said simulation, 
calculating individual sound signals to be emitted by at 
least some of the transducers of said array of 
transducers in order to obtain said target soundwave 
field; 

• and a stage of actually positioning the array of 
transducers on the brain barrier, this stage comprising 
the following steps: 
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2a) initially positioning the array (17) of 
transducers on the bone barrier (14) approximately in the 
specific position; 

2b) using at least some of the transducers of 
the array of transducers to perform echography to locate 
the position of said array (17) of transducers relative 
to the bone barrier (14) ; and 

2c) refining the position of the array (17) of 
transducers relative to said bone barrier as a function 
of the location determined in step 2b) . 

6. A method according to any preceding claim, in which 
the soundwaves are at frequencies lying in the range 

0 - 5 MHz to 3 MHz . 

7. A method according to any preceding claim, in which 
step le) is followed by a step If) during which emission 
of sound signals by the array (17) of transducers is 
simulated, said signals being determined from said 
individual sound signals and serving to obtain a desired 
soundwave field, propagation of soundwaves generated by 
said emission is simulated, and it is verified that said 
propagation satisfies certain predefined criteria. 

8. A method according to any preceding claim, in which: 

■ during step Id) , propagation of soundwaves from at 
least one point (18) in the substantially homogeneous 
medium towards at least some transducers of the array of 
transducers is simulated, and received simulated sound 
signals Ri(t) reaching the locations of said transducers 
± of the array of transducers are determined, where i is 
an integer in the range 1 to n, and n is the number of 
transducers in the array of transducers; and 

• during step le) , the individual sound signals 
Ei(t) for emission by each transducer i under 
consideration are determined as being proportional to a 
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time reversal Ri(-t) of said received simulated sound 
signals Ri (t) as previously determined in step Id). 

9. A method according to claim 8, in which, during step 
le) , the sound signals Ei (t) to be emitted are determined 
by the formula: 

Ei (t) = Gi .Ri (-t) 
where Gi is a gain factor that differs from one 
transducer i to another, for compensating dissipation in 
the bone barrier. 

10. A method according to claim 9, in which the gain 
factors Gi corresponding to at least some of the 
transducers are respectively inversely proportional to 
the square of an amplitude of the corresponding received 
simulated sound signals Ri (t) . 

11. A method according to claim 8, in which: 

- during step Id) , the simulation is performed by 
using a virtual three-dimensional map of soundwave 
absorption, having absorption coefficients -x at each 
point of the bone barrier that are opposite to the real 
absorption coefficients x determined during step lb) ; and 

• during step le) , the individual sound signals 
Ei(t) to be emitted are determined as being equal to said 
time reversal Ri(-t). 

12. A method according to any one of claims 1 to 7, in 
which the array (17) of transducers is included in a 
fluid- filled tank (13) having at least one flexible wall 
(12) for pressing against the bone barrier (14) , the 
specific positions of the transducers as determined 
during step Ic) not being in contact with the bone 
barrier, and in which: 

during step Id) : 

• propagation of soundwaves from at least one point 
(18) of the substantially uniform medium (19) towards at 



least some transducers of the array (17) of transducers 
is determined, and received simulated sound signals Ri(t) 
reaching the locations of said transducers of the array 
of transducers are determined, where ± is an integer in 
the range 1 to n, and where n is the number of 
transducers in the array of transducers; 

- then emission by each transducer i^ of a sound 
signal Ri(-t) corresponding to a time reversal of the 
signal Ri (t) is simulated, and propagation in said fluid 
to a virtual transducer i_ situated in contact with the 
bone barrier (14) in correspondence with the transducer ±_ 
is simulated, and received simulated sound signals R»i(t) 
reaching the location of said virtual transducer _i are 
determined; 

• then emission by each virtual transducer ±^ of an 
acoustic signal G'i.R'i(-t) is simulated where R'i(-t) is 
a time reversal of the signal R*i(t) and where G'i is a 
reversal coefficient proportional to the square of an 
amplitude of the signal R'i(t), at least for some of the 
virtual transducers i_; 

• then propagation in said fluid to the transducer i 
is simulated and received simulated sound signals R"i(t) 
reaching the location of said transducer ± are 
determined ; 

and during step le) , the individual sound signals Ei(t) 
to be emitted are determined as being equal to a time 
reversal R"i(-t) of said received simulated sound signals 
R"i (t) . 

13. A method according to any one of claims 1 to 7, in 
which during step Id) , the emission of a soundwave pulse 
by at least some of the transducers ± of the array (17) 
of transducers is simulated and propagation of soundwaves 
from each transducer _i in consideration towards a 
plurality of reference points r situated in the 
substantially uniform medium (19) is simulated, where i 
is an index in the range 1 to n designating a transducer 
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of the array, and where n is a non-zero natural integer 
designating the number of transducers, r being an integer 
in the range 1 to m, where m is a non-zero natural 
integer designating the number of reference points, and 
5 simulated impulse responses hri(t) reaching each of said 
reference points r of the substantially homogeneous 
medium are determined, step le) comprising the following 
substeps : 

lei) determining a number g of frequency components 
10 for each of the simulated impulse responses, having 

respective frequencies cok, where k is an index lying in 
the range 1 to p and designating a frequency component; 

le2) determining p transfer matrices H(cok) = 
[Hri (cok) ] , _i lying in the range 1 to n, and r going from 
15 1 to m, where Hri (cok) is the value at the frequency cok of 
the Fourier transform of the impulse response Hri(t); and 

le3) for each reference point r, n components 
Ei(cok,r) are determined such that: 

F(<ok,r) = H(ci)k) .E(cok,r) 
20 where E(cok,r) = [Ei(cok,r)] is a vector having n 

components Ei(a)k,r), F(a)k,r) = [F^((ok,r)] is a vector 
having m components F^(cok,r) where I varies in the range 1 
to m, these m components F^(cok,r) corresponding to 
generating said predetermined target soundwave field at 
2 5 the frequency cok at the points r. 

14. A method according to claim 13, in which, during 
substep le3) , g matrices H'l (cok) are calculated at least 
by inverting the transfer matrices H(cok), and for each 
30 reference point r of the substantially homogeneous 
medium, the vector E(cok,r) is calculated using the 
formula: 

E(cok,r) = H-Mcok) .F(cok,r) 

35 15. A method according to any preceding claim, in which, 
during step Id), the impulse responses hir(t) between a 
plurality of reference points r of the substantially 



homogeneous medium and at least some of the transducers i 
of the array (17) of transducers are determined, where i 
is an index in the range 1 to n which designates a 
transducer, and where n is a non-zero natural integer 
which designates the number of transducers, r being an 
integer lying in the range 1 to m, m being a non-zero 
natural integer designating the number of reference 
points, and during step le) , it is also determined how to 
focus at least a portion of the array (17) of transducers 
in reception on each reference point r in order to make 
an echographic image. 

16. A method of medical imaging by echography, the method 
comprising a non- invasive method of obtaining a target 
soundwave field according to any preceding claim, and an 
imaging stage during which at least one echographic image 
of the substantially homogeneous medium (19) is made 
using at least a portion of the array (17) of 
transducers, using the individual sound signals as 
determined during the training stage. 

17. Apparatus specifically designed for implementing a 
method according to any preceding claim, the apparatus 
comprising at least: 

• an array (17) of transducers adapted to be 
positioned outside a bone barrier (14) masking a 
substantially homogeneous medium (19) ; 

• mapper means (7) for determining three-dimensional 
maps of at least density, soundwave speed, and soundwave 
absorption in the bone barrier, on the basis of a three- 
dimensional image of said bone barrier made by X-rays and 
giving the porosity of said bone barrier at each point; 

• simulator means (7) for simulating at least one 
propagation of soundwaves between at least one point (18) 
of the substantially homogeneous medium and at least some 
of the transducers of the array (17) of transducers on 
the basis of a mathematical model of propagation and on 



the basis of said three-dimensional maps of density, of 
soundwave speed, and of soundwave absorption, and as a 
function of a specific position for the array (17) of 
transducers relative to the bone barrier (14) ; 

• calculator means (7) for responding to said 
simulation to calculate individual sound signals to be 
emitted by at least some of the transducers of said array 
(17) of transducers in order to obtain a target soundwave 
field in the substantially homogeneous medium; 

• locator means (7) for using at least some of the 
transducers of the array (17) of transducers to locate 
the position of said array (17) of transducers and 
relative to the bone barrier (14) by echography; and 

• position-refiner means (7, 4) for refining an 
initial position of the array (17) of transducers 
relative to the bone barrier (14) as a function of the 
position of the array of transducers as located relative 
to the bone barrier in such a manner that the position of 
the array of transducers relative to the bone barrier 
corresponds to the specific position. 

18. Apparatus according to claim 17, in which the means 
for locating the position of said array (17) of 
transducers relative to the bone barrier (14) are adapted 
to determine an outside shape of at least a portion of 
the bone barrier (14) by echography, using at least a 
portion of said array (17) of transducers, by comparing 
said outside shape with said three-dimensional image of 
said bone barrier. 

19. Apparatus according to claim 17 or claim 18, 
comprising a locating device (8) provided with securing 
means (16) adapted to secure said locating device rigidly 
on the bone barrier (14) , said locating device being 
adapted to absorb X-rays and being visible in the three- 
dimensional image of said bone barrier, and the means (7) 
for locating the position of said array (17) of 



transducers relative to the bone barrier (14) are adapted 
to locate the position of the array (17) of transducers 
relative to the locating device (8) by echography. 

20. Apparatus according to any one of claims 17 to 19, in 
which the array (17) of transducers is included in a 
fluid-filled tank (13) having at least one flexible wall 
(12) for pressing against the bone barrier (14) . 

21. Apparatus according to any one of claims 17 to 20, in 
which the array (17) of transducers comprises both a sub- 
array (17b) for imaging and a sub-array (17a) for 
hyperthermia treatment, these two sub-arrays comprising 
transducers of respective different types. 

22- Apparatus according to any one of claims 17 to 21, in 
which the array (17) of transducers is adapted to emit 
soundwaves at frequencies lying in the range 0.5 MHz to 
3 MHz- 

23. Apparatus according to any one of claims 17 to 22, 
further comprising simulator means (7) adapted: 

• to simulate at least the array (17) of transducers 
emitting sound signals determined from said individual 
sound signals and enabling a desired soundwave field to 
be obtained; 

• to simulate propagation of soundwaves generated by 
said emission; and 

• to verify that said propagation satisfies certain 
predefined criteria. 

24. Apparatus according to any one of claims 17 to 23, in 
which : 

the simulator means (7) are adapted to simulate 
soundwave propagation from at least one point (18) of the 
substantially homogeneous medium towards at least some of 
the transducers of the array of transducers to determine 



received simulated sound signals Ri(t) reaching the 
locations of said transducers i of the array of 
transducers, i being an integer in the range 1 to n, and 
n being the number of transducers in the array of 
transducers; and 

the calculator means (7) are adapted to determine 
the individual sound signals Ei (t) to be emitted by each 
transducer ± under consideration as being proportional to 
a time reversal Ri(-t) of said received simulated sound 
signals Ri (t) as previously determined in step Id). 

25. Apparatus according to claim 24, in which the 
calculator means (7) are adapted to determine the sound 
signals Ei(t) to be emitted using the formula: 

Ei (t) = Gi .Ri (-t) 
where Gi is a gain factor that differs from one 
transducer i to another, in order to compensate for 
dissipation in the bone barrier. 

26. Apparatus according to claim 25, in which the gain 
factors Gi corresponding to at least some of the 
transducers are respectively inversely proportional to 
the square of an amplitude of the corresponding received 
simulated sound signals Ri (t) . 

27. Apparatus according to claim 24, in which: 

the simulator means (7) are adapted to perform said 
simulation of soundwave propagation by using a virtual 
three-dimensional map of soundwave absorption having at 
each point of the bone barrier absorption coefficients of 
-T opposite to the real absorption coefficient x as 
determined during step lb); and 

the calculator means (7) are adapted to determine 
the individual sound signals Ei (t) to be emitted as being 
equal to said time reversal Ri(-t). 



28. Apparatus according to any one of claims 17 to 23, in 
which the array (17) of transducers is included in a 
fluid-filled tank (13) having at least one flexible wall 
(12) for pressing against the bone barrier (14) , the 
locations provided for the transducers taken into account 
by the simulator means (7) not being in contact with the 
bone barrier, 

in which apparatus the simulator means (7) are 
adapted : 

• to simulate propagation of soundwaves from at 
least one point (18) of the substantially homogeneous 
medium (19) towards at least some of the transducers of 
the array (17) of transducers and to determine the 
received simulated sound signals Ri(t) reaching the 
locations of said transducers i of the array of 
transducers, i being an integer in the range 1 to n, and 
n being the number of transducers in the array of 
transducers ; 

• to simulate emission by each transducer _i of a 
soundwave Ri(-t) corresponding to a time reversal of the 
signal Ri (t) , and to simulate propagation in said fluid 
to a virtual transducer ± situated in contact with the 
bone barrier (14) in correspondence with the transducer 
i, and to determine the received simulated sound signals 
R'i(t) reaching the location of said virtual transducer 

i ; 

• to simulate emission by each virtual transducer i 
of a sound signal G'i.R'i(-t) where R'i(-t) is a time 
reversal of the signal R'i(t), and where G'i is a 
reversal coefficient proportional to the square of an 
amplitude of the signal R'i(t), at least for some of the 
virtual transducers and • 

• to simulate propagation in said fluid to the 
transducer i, and to determine the received simulated 
sound signals R"i(t) reaching the location of said 
transducer i; 



and in which the calculator means (7) are adapted to 
determine the individual sound signals Ei(t) to be 
emitted as being equal to a time reversal R"i(-t) of said 
received simulated sound signals R"i(t)- 

29. Apparatus according to any one of claims 17 to 23, in 
which the simulator means (7) are adapted: 

• to simulate emission of a soundwave pulse by at 
least some of the transducers i of the array (17) of 
transducers, and propagation of soundwaves from each 
transducer ± in question towards a plurality of reference 
points r situated in the substantially homogeneous medium 
(19) , ± being an index in the range 1 to n designating a 
transducer of the array, and n being a non-zero natural 
integer designating the number of transducers, r being an 
integer. in the range 1 to m, and m being a non-zero 
natural integer designating the number of reference 
points ; 

• to determine simulated impulse responses hri(t) 
reaching each of said reference points r of the 
substantially homogeneous medium; 

and in which the calculator means (7) are adapted: 

• to determine a number p of frequency components 
for each of said simulated impulse responses at 
respective frequencies cok, k being an index in the range 
1 to £ and designating a frequency component; 

• to determine £ transfer matrices H(ook) = 

[Hri ((Ok) ] , i lying in the range 1 to n and r lying in the 
range 1 to m, where Hri (cok) is the value at the frequency 
cok of the Fourier transform of the impulse response 
Hri (t ) ; and 

• to determine, for each reference point r, n 
components Ei(cok,r) such that: 

F(cok,r) = H(cok) .E(cok,r) 
where E(cok,r) = [Ei((ok,r)] is a vector having n 
components Ei(cok,r), F(cok,r) = [F-^(cok,r)] is a vector of 
m components F^(cDk,r) where £ lies in the range 1 to m. 
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these m components F-^((Dk,r) corresponding to generating 
said predetermined target soundwave field at the 
frequency cak at the points r. 

30. Apparatus according to claim 29, in which the 
calculator means (7) are adapted: 

" to calculate p matrices H-=^ (ok) at least by 
inverting the transfer matrices of H (cok) ; and 

• for each reference point r of the substantially 
homogeneous medium, to calculate the vector E(cok,r) by 
means of the formula: 

E(cok,r) = H-i(cok) .F(a)k,r) 

31. Apparatus according to any one of claims 17 to 30, in 
which the simulator means (7) are adapted: 

• to determine the impulse responses hir(t) between 
a plurality of reference points r of the substantially 
homogeneous medium and at least some of the transducers i 
of the array (17) of transducers, ± being an index lying 
in the range 1 to n designating a transducer, and n being 
a non-zero natural integer designating the number of 
transducers, r being an integer lying in the range 1 to 
m, and m being a non-zero natural integer designating the 
number of reference points; 

and in which the calculator means (7) are adapted to 
determine how to focus at least a portion of the array 
(17) of transducers in reception on each reference point 
r in order to make an echographic image. 

32. Apparatus according to any one of claims 17 to 31, 
further comprising imaging means (7, 7b) adapted to make 
at least one echographic image of the substantially 
homogeneous medium (19) using at least a portion of the 
array (17) of transducers, using the individual sound 
signals determined by the calculator means (7) . 



